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Abstract 

The question of constructing the finite genus quasiperiodic solutions for the Ablo- 
witz-Ladik hierarchy (ALH) is considered by establishing relations between the ALH 
and the Fay's identity for the ^-functions. It is shown that using a limiting procedure 
one can derive from the latter an infinite number of differential identities which can 
be arranged as an infinite set of differential-difference equations coinciding with the 
equations of the ALH, and that the original Fay's identity can be rewritten in a form 
similar to the functional equations representing the ALH which have been derived in 
the previous works of the author. This provides an algorithm for obtaining some class 
of quasiperiodic solutions for the ALH, which can be viewed as an alternative to the 
inverse scattering transform or the algebro- geometrical approach. 

1 Introduction. 

The problem of constructing the quasiperiodic solutions (QPS) is one of the most challeng- 
ing problems of the theory of integrable systems, and many mathematicians and physicists 
spent much efforts to obtain the QPS for almost all equations that are known to be in- 
tegrable. The Ablowitz-Ladik hierarchy (ALH), which has been introduced in is not 
an exception. So, e.g, one should mention the works by N.N.Bogolyubov (Jr) et al 0. § 
and S.Ahmad, A.R.Chowdhury [|, H devoted to the discrete nonlinear Schrodinger equation 
(DNLSE) and the discrete modified Korteveg-de Vries equation (DMKdV), which are the 
best studied equations of the ALH. There authors were studying this problem in the frame- 
work of the inverse scattering transform (1ST). Another, the so-called algebraic-geometrical, 
approach has been used by Miller, Ercolani, Krichever and Levermore, who considered in 
H the complex version of the DNLSE and obtained the Baker-Akhiezer function and QPS 
corresponding to finite genus Riemann surfaces. This work provides almost exhaustive solu- 
tion of the problem of the finite genus QPS, but its results need some further simplification 
to be useful for practical purposes, especially if one wants to extend them to the higher 
equations of the ALH and in this work I will try to avoid the algebro-geometrical language, 
and will use some more direct (and simpler) strategy. As it has been established in [[J, 
each finite genus QPS of the DNLSE can be presented as a quotient of the ^-functions of 
some arguments multiplied by an exponent of some phase, all of them being some linear 
functions of the coordinates (the same is true in the cases of the DMKdV as well as all other 
equations of the hierarchy). Thus, since we know the structure of the solutions, all we have 
to do to derive them is to calculate some number of constant parameters. So, it desirable to 
develop some method which will enable us to obtain these constants (and hence solutions) 



straightforwardly, not using technique (sometimes rather complicated) of the theory of the 
functions and differentials on hyperelliptic Riemann surfaces. It turns out that this can be 
done. Moreover, this can be done not only for the DNLSE or DMKdV but, in principle, for 
all equations of the hierarchy simultaneously. Namely this is the main question of this paper. 
The key moment is that the ^-functions of the finite genus Riemann surfaces (of which the 
finite genus QPS are built-up) satisfy some algebraic relation, the so-called Fay's trisecant 
formula 0, 0], which can be used to obtain an infinite number of differential identities, which, 
as will be shown below, are closely related to the ALH, and can be used to obtain the QPS 
we are looking for. Such approach also demonstrates some new, to my knowledge, feature of 
the ALH (and namely this was one of the main motives to write this paper): the equations 
of the ALH naturally appear when flows over Riemann surfaces are considered (I will return 
to this question below). 

The plan of the paper is as follows. After presenting some basic facts on the ALH (section 
H) I will discuss the Fay's formula and its differential consequences (section |). These results 
will be used to obtain the finite genus QPS for the ALH (section f|). 



2 Ablowitz-Ladik hierarchy 

The ALH is an infinite set of integrable differential-difference equations, which has been 
introduced in All equations of the ALH can be presented as the compatibility condition 
for the linear system 



# n+1 = U n V n (2.1) 

^vD„ = V ^ n , j = ±l,±2,... (2.2) 

where *ff n is a 2-column, U n and V n are 2x2 matrices with U n being given by 

U n = U n {\)=(^ M (2.3) 



Qn A~ 

(here A is the auxiliary (spectral) constant parameter) and the matrices Vn are polynomials 
in A, A -1 . The ALH in a natural way can be splitted into two subsystems (subhierarchies) . 
One of them corresponds to the case when Vn , j = 1, 2, . . . are jth order polynomials in 
A -1 ('positive' subhierarchy) . Its simplest equations are 

— = -iPnQn+i (2-4) 

-7^ = iPnTn-l (2.5) 

where 

p n = 1 - q n r n (2.6) 

The 'negative' subhierarchy is build up of the ^-matrices being polynomials in A. Its simplest 
equations are 
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dq n 
dzi 
dr n 
dz\ 

(I use the notation Zj = Z-j, j = 1, 2, . . .). 



-ip n q n -x 



(2.7) 
(2.8) 



It has been shown in |T(| that the ALH can be presented in the form of the functional- 
difference equations: 



q n (z,z) - q n (z - i{£],z) = £ 
r n (z,z)-r n (z + i[ti],z) = £ 
for the 'positive' subhierarchy, and 



1 - q n (z,z)r n (z - i[£],z) 
1 - q n (z + i[£],z)r n (z,z) 



q n +i{z,z) 
r n -i(z,z) 



(2.9) 
(2.10) 



q n (z,z) -q n (z,z-i[£ *]) 
r n (z,z) - r n (z,z + 



-i 



l-q n (z,z)r n (z,z-i[£ ]) q n -x(z,z) (2.11) 



C 1 1 -q n (z,z + i[C l ])r n (ZiZ) r n+1 (z,z) (2.12) 



for the 'negative' one. Here the designations 



and 



f(z,2±[r 1 }) 



f(z, z) — f [Zi, z 2 , z 3 , . . . z±, z 2 , z 3 , . . .) 



/ ( Zl ±£,z 2 ± e/2, z 3 ± £ 3 /3, ...2^22, z 3 , . . .) 

f (zi, z 2 , z 3 ,...z l ± r \ z 2 ± r 2 /2, zs ± r 3 /s, . . .) 



(2.13) 



(2.14) 
(2.15) 



are used. Expanding equations Q2.9Q , (|2.10| ) in power series in £ one can obtain all equations 
of the 'positive' subhierarchy. Analogously, expanding equations ( 2.11] ), ( |2.12j ) in power 
series in £ _1 one can obtain all equations of the 'negative' one. 

In what follows I will use also the tau-functions of the ALH, a n , p n and r„, which are 
defined by 



On 

q n = — , 



Pn 



(2.16) 



and 



T~n-lT n +l =T%- O n p n (2-17) 

The functional representation of the ALH in terms of the tau-functions can be written as 



r n (z) a n (z + «[£]) - a n (z) r n (z + 

Pn{z) T n (z + - T n (z) p n (z + 

r n (z) T n (z + - p n (z) a n (z + 



^r n ^{z)a n+1 {z + m) (2.18) 
Cp n ^{z)r n+1 {z + i[^)) (2.19) 
T n ^(z)T n+1 (z + i^)) (2.20) 
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(where dependence on z/s is omitted) and 



r n (z) a n (z + zfr 1 ]) - <T n (z) r n (z + ^r 1 ]) = r 1 r n+1 (z) o n _ x (z + z^" 1 ]) (2.21) 

Pn (z) rn^+ifr 1 ]) -r„(f) Pn^+itr 1 ]) = r'pn+i^r^i^+^r 1 ]) (2.22) 

r„(z) r n (z + ^r 1 ]) - Pn(z) a n (z + i^ 1 ]) = r n+1 (z) r n ^(z + z^ 1 ]) (2.23) 

(where dependence on Zj's is omitted). 

The key idea of the present work is to establish the relation between these equations and 
the famous Fay's identity for the ^-functions which can be used to derive the finite-gap QPS 
of the ALH. 



3 Fay's identity. 

In this paper we will deal with the compact Riemann surface X of the genus g corresponding 
to the hyperelliptic curve 

s 2 =V 2g+2 (0 (3.1) 

where V 2g+2 (^) is a polynomial without multiple roots of degree 2g + 2. In the framework 
of the 1ST such curves appear in the analysis of the scattering problem Q2.1Q . For example, 
in the case of the periodic conditions 

Qn+g+l = Qnj r n+g+l = r n (3-2) 

the polynomial 7 3 2g+2(0 is defined by 

V 2g+2 (\ 2 ) = A^+D {[tr T„(A)] 2 - 4detT n (A)} (3.3) 
where T n (A) is the transfer matrix of the scattering problem ( |2 . 1| ) . 

T n (X) = U n+g (X)...U n (X) (3.4) 

(it can be straightforwardly shown that the right-hand side of ( |3.3|) under the restriction 
( p.2| ) does not depend on the index n). Topologically, X is a sphere with g handles. One 
can choose a set of 2g closed contours (cycles) {a,, &i}i=i,... )g with the intersection indices 

cii o ctj = bi o bj = 0, di o bj = 5ij i, j = 1, . . . , g (3.5) 
and find g independent holomorphic differentials, say ones given locally by 

Uk = k = l,...,g (3.6) 

which can be used to construct the canonical basis of the holomorphic 1-forms 

9 

^fe = ^ C k>l ut (3.7) 
i=i 

where u^s satisfy the normalization conditions 
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u k = 5 ik (3.8) 
Then, the matrix of the 6-periods, 

Vt ik = (h uj k (3.9) 



b, 



determines the so-called period lattice, = {m + Qn, m, n G Z 5 }, the Jacobian of this 
surface Jac(X) = C 9 /L n (2g torus) and the Abel mapping X — > Jac(X), 

P ^ I u> (3.10) 

Jpo 

where us is the ^-vector of the 1-forms, us = (ui, . . . , uj g ) T and Pq is some fixed point of X. 
A central object of the theory of the compact Riemann surfaces is the ^-function, 9(C) — 

e(c,n), 

9(C) = ex P{^ nfin +2mnC} (3.11) 
which is a quasiperiodic function on C 9 



9(C + n) = 9(C) (3.12) 
9(C + Qn) = exp{-mnQn -2ninC}9(C) (3.13) 

for n G Z 9 . 

To simplify the following formulae I will use the designations 

o A B (C) =6 (t + J* u ) ^ 3 - 14 ) 

and 

9 A B =0[*[(j A (3-15) 
Here 9 [c] (C) is the so-called ^-function with characteristics, 

9 [c] (C) = exp {iri atta + 2m a (C + b)} 9 (C + tta + b) , c=(a,b) (3.16) 
and 8 = (6', 6") G |Z 2fl /Z 29 is a non-singular odd characteristics, 

9 [6] (0) = 0, grad c 9 [6] (0) ^ (3.17) 
Function fl3.15| ) is closely related to the prime form [[7], || . 



On 

E(P, Q) = Q - = (3.18 

VW)v^(Q) 



where \ is given by 
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^) = E(|-^[^])(0)^(P) (3.19) 

The prime form E(P, Q) is skew-symmetric, E(P, Q) = —E(Q, P), has first-order zero along 
the diagonal P = Q and is otherwise non-zero. Analogously, 

0q = -e Q p , e p p =o (3.20) 

One of the most interesting results of the theory of the ^-functions is the following identity 
for the ^-functions associated with the finite-genus Riemann surfaces, the Fay's identity: 

e p P \ K\ 0£(C) *2(C) - K\ K\ e p P %) K%) = K\ 6% d(C) 0%%{C) (3.2i) 

(here Pi,...,P± are arbitrary points of X) and namely this formula will be basis of the 
following consideration. 



4 Quasiperiodic solutions. 

It is already known that in the quasiperiodic case the tau-functions of the ALH are (up to 
some simple factors) the ^-functions of different arguments and I am going now to present 
the Fay's identity and some of its corollaries in the form similar to fl2.18| )-( P?.20| ) and ( [2.21 )- 



( p.23| ) which will enable to obtain the finite-gap solutions of these functional equations, i.e. 
to obtain the finite-gap solutions of the ALH. 

Hereafter I will use the letters A, B, C, D for the points of the Riemann surface which 
correspond to the points and oo of the complex plane, 

A = oo + , D = oo_, P = 0_, C = 0+ (4.1) 

Since A,D and B, C are poles and zeroes of the meromorphic function 7r(P), which is a 
projection of X onto the extended complex plane P 1 sending a point P = (s, £) into £, they 
satisfy, according to the Abel's theorem, the condition f B ^ u) G Lq. The integration pathes 
in ( |3.10|) can be chosen in such a way that 

AD 

u> = (4.2) 

BC 



(here zero stands for from C 9 , not from Jac(A)) and in what follows I will accept (|4.2p as 
true. 

Now I am going to use ( |3.21| ) thinking of three points from (Pj_, P 2 , P3, P4) as constant 
(I will choose them from the set (A, B, C, D)) and the fourth one (I will denote it by P) as 
variable. Setting (Pi, P 2 , P3, P4) = (A, B,C, P) one can rewrite ( |3.21| ) as 

o B a K 0(C) P D (C) + H o p c e A c {C) e p B (C) = h % ^(0 e p c (C) (4.3) 

This formula is the quasiperiodic analog of ( |2.18| ). Shifting the arguments of the ^-functions, 
C — ► C + J A u>, one can obtain the equation which will be transformed below to ( 2.19| ): 



o B c o p A em em + H K e(C) e c /m = e A c e p em em (4.4) 
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At last, replacing in (|3~2lD (P 1; P 2 , P 3 , P 4 ) with (A, P, C, D) using flOj) and making the shift 
£ — > C + X4 u one can write the identity 

o A J p c 9(C) o p B (C) -e A J p A m) m) =o A c e p D e A B {t) e p A (C) (4.5) 

which is a quasiperiodic analogue of ( |2.20| ). 

Our first goal is to present equations ( |4.3|) -( |475|) in the bilinear form. To this end I will 
first shift the arguments of the ^-functions: £ — > £ n , 

C n = C + n [ u (4.6) 

J A 

Next, I will introduce the functions a n (P), p n (P) and r n (P), 



r n (P) = a n (P)9 p B (U (4.7) 

a n (P) = [3 n (P)6 p D (U (4.8) 

Pn (P) = 7 n(P)C(0 (4-9) 

It is not difficult to verify that if one chooses the functions a n , /3 n , j n as follows 



a n {P) = a*// 1 I 2 exp{rup DC (P)} (4.10) 
/3 n (P) = q*e n exp{<p AC (P)}a n (P) (4.11) 
7 „(P) = r*e- n exp{-<^ (P)}a n (P) (4.12) 



where the functions (p QR are defined in the vicinity of the point B by 



exp {cp QR (P) - cp QR (B)} = Q -fS- (4.13) 



the constant \i is given by 

H = \^- (4.14) 
and a*, g#, r* and e are arbitrary constants satisfying 



2 

9 A 



2 



u D U C 

then (FO)-(PO) can be rewritten in terms of the functions a n (P), p n {P) and T n (P) as 



r n (B) a n (P) - a n {B) r n (P) = K(P) r n ^(B) a n+1 (P) (4.16) 
p n (P)r n (P)-r n (P)p n (P) = K(P) Pn -i(B) r n+1 (P) (4.17) 
x„(P)x„(P)-p n (P)a n (P) = r n _a(P) r n+1 (P) (4.18) 



where 
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1 OA OP OP 

KIP) = - (4.19) 

U C U A u D 

Thus we have presented the Fay' identities in the bilinear form similar to ( |2.18 )- (|2.20| ). 
What I have to do now is to introduce a z-dependence in such a way that a shift over the 
Riemann surface from the point B to a point P (which correspond to the points and £ of 
the complex plane) can be taken into account by the simultaneous shifts z m — > z m + z£ m /m: 



f n (B) = f n (z) = f n (z h z 2 ,z 3 ,...), (4.20) 

fn{P) = fn(z + t[Q) = fn(Zl+t^Z 2 +te/2,Z 3 + te/3,...) (4.21) 

(I hope that the usage of the same letters for functional dependence on both the point of 
the Riemann surface and the ALH variables z m will not lead to confusion). In other words 
I want to introduce such functions C(zi, z 2 , ■■■) and f QR (zx, z 2 , ■■■) that 

r-P 



C(z + i[e])-C(*)= / " (4.22) 
Jb 



and 



¥>qr(z + ~ Vqr(z) = Vqr{P) ~ <Pqr(B) (4.23) 

This can be done as follows. In the neighborhood of the point B (which is an preimage of the 
point £ = of the complex plane) the components of the integral in (|4.22 ) can be presented 
in terms of the ^-coordinate as 

Jb J-f Jo \/V 2g+2 {x) 

where the sign of the square root is fixed by vT = I. Hence, taking £ to be a linear function 
of the coordinates z m , 



oo 



c = cc*) = Yl s {m) z m ( 4 - 25 ) 

m=l 

one can conclude that to satisfy ( f4.22| ) the vectors C should be defined as the coefficients 
of the series 

OO j 

^ c M r = _- e W (0 (4.26) 

m=l 



(here W is the vector with the components Wk). Using ( 4.24 ) one can rewrite ( 4.26|) as 



oo 9 ti 

E<Tf = /v k (4.27) 

m=l 1=1 V ' 29+2 IS I 

(the right-hand side of this equation should be understood as a power series in £). 

In a similar way one can tackle equation ( |4.23|) and to derive the following result: ip QR (z) 
is the linear function, 



S 



m=l 



with the coefficients being defined by 

(the right-hand side is again a power series in £). 

Thus one can write the following expressions for the tau-functions: 



(4.28) 



(4.29) 



T n (z) = a*/i" 2/2 exp{ mp DC (z)} 9 (C n (z)) 



<Jn(Z) 



Pn{Z) 



C 



(4.30) 
(4.31) 



q*e n fi n 11 exp { n tp DC (z) + <p AC (z)} 9 ( C„0) ~ J " 
ne- n fx n2 / 2 ex P {n^ DC (z)-^ AC (z)} dU n {z)+ J^e^ (4.32) 



At last, we have to rewrite the function K(P) from the right hand side of ( |4.16|) - (|4.18|) . 
This is the first time, since the Fay's identity has been written down, that we need some 
facts from the theory of the Riemann surfaces - till now all was done by simple algebra. 
Consider the function 



f(P) 



B U C 

9 P 9 P 



(4.33) 



This is a single- valued, due to the condition ( |4.2|) , function which possesses zeroes at the 
points B, C and poles at A, D. Remembering that B, C correspond to £ = 0, and A, D - 
to £ = oo, one can easily obtain one function with the same divisor, B + C — A — D, namely, 
the projection ir(P) discussed above (see paragraph before ((4.2|)). The quotient ir(P)/ f(P) 
has no poles (and zeroes as well) on X, hence it is a constant 



Thus, if we take 



f(P) = C£ for P=(s,0 



e = C 6 -X 
B c 



(4.34) 



(4.35) 



then 



K{P) = £ 



(4.36) 



and relations ( |4.16 )-(4.18) become ( ^.18[ )-( p.2(J| ), or in other words, the functions defined by 
(CT)-(ra) solve equations (|2~T^)-(|2T20|). 

Till now we were operating in a neighborhood of the point B and obtained solutions 
of equations ( |2.18|) - (|2.20| ), and hence of ( |2.9| )- (|2.10|) , i.e. solved the 'positive' part of the 
ALH. To take into account the 'negative' equations (pHD - (pl2|) , or (|2T2l|) - (|2~23D , one can 
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proceed in the similar way, but this time considering flows near another distinguished point, 
D, which is an preimage of the point £ = oo. It can be shown that functions r n , a n , p n 
given by ( }4.30| )-( f4.32| ) will solve fl2.21| )-( p.23j ) provided we introduce the z-dependence by 
replacing 



C(z) - C(z,z) (4.37) 

¥dc{z) -> ip DC (z) + if BA (z) (4.38) 

Vac{z) -> <Pao{z) + <Pca(z) (4-39) 

(the overbar does not mean the complex conjugation!) where 



C(z,z + iiC 1 })-C(z,z) = W(C 1 )= f u (4-40) 

JD 



and 



<Pq*(* + i[C 1 ]) ~ <P Q n{z) =^ e f- 9 -l (4.41) 

Thus, now we have all necessary to formulate the main result of this paper. The finite genus 
solutions of the ALH can be presented as 

q n (z,z) = g , e »exp {*,(*,*)} H %;Pt r f~y ) (4-42) 

9 {Q(z, z) + nil) 

m 0(C(z,z) +nU + V) lt in . 



'B pC 

U = I w, V = / u (4.44) 



where 

J7 = 

J A J A 

The functions z) and (^(2;, 2) are given by 

00 

C(z, z) = (^ (m) z ™ + z™) + constant (4.45) 

m=l 

00 

^(z, z) = {V (m) z m + ¥ {m) z m ) + constant (4.46) 

m=l 

where the constants C^ m \ C an d <^ m ) are defined as coefficients of the series 



m=l 1=1 V ' 29+2 Is J 

E5 m) r m = -igc?^-,^— (4.48) 
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with 



TWO = £ 2s+2 (£)7V 2 (i/£) 



(4.49) 



and 



E ^ (m) £ r 

m=l 

oo 

E ^ (m) r r 



d. fl[3](C/- y + W(Q) 

H de e[(5](c/+w(e)) 

d 0[*](tf + W(i/o) 



m=l 



*d£ 6[5)(U + V + W{l/0) 



(4.50) 
(4.51) 



The constant e is given by (|4.35|) and Q* , are arbitrary constants related by ([4.151) . 
The 'real' tau-function r„ can be written as 

r n (z, z) = a,/i" 2/2 exp { n ^{z, z)} 6 (C(z, z) + nU) (4.52) 
where the constant /i is given by (|4.14|) , 



ijj( Z} z) = E (V' (m) z m + ^ (m) ^n) + constant 



m=l 



and 



(4.53) 



E^ (m) £ r 

m=l 

oo 

E^ V- ,(m ' ) 



d_ . g[^](c/- v + w(Q) 

? de 0[«5](-v+w(o) 



m=l 



e [6] (v + wwq) 

d£ *"#[<5] (|7 + V + W(l/0) 



z ^ — — In 



(4.54) 
(4.55) 



5 Discussion. 

In this paper we have obtained the finite genus solutions for the ALH. The results can be 
used to derive also the finite genus solutions for other integrable hierarchies, which can be 



embedded' into the ALH (see [11, 12, 10, 13]. So, for example, the functions 



where 



e (CO*, 2) " v) 

Q = — = Q*exp{ip(z,z) +^(z,z)} — — 

r u {(,{z, z)) 

eU(z,z) + v 

R = =R m exp{-i(>(z,z)-<p(z,z)} — 77-7 — 

T V{^[Z,Z)) 



V = V - u 

and constants Q*, R* are related by 



c 



UJ. 



B 



(5.1) 
(5.2) 



(5.3) 
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solve the nonlinear Scrodinger equation 



QA OA 
9 A 9 B 



(5.4) 



id 2 Q + d n Q + 2Q 2 R = 
-id 2 R + d n R + 2QR 2 = 



(5.5) 
(5.6) 



where d m = d/dz m , as well as all higher equations of the hierarchy (see fllOl). 
The quantities 

Tn+lTn-l = {CA^ x) + (n - 1)U) 9 (C T (x, x) + (n + l)U) 

where 

Ct{ x -> x ) = C^ x + C^ x + constant 
solve the 2D Tocla lattice equation 

d 2 

-\np n = p„_i - 2p n + p n+ i 



dxdx 



(5.7) 



(5.8) 



(5.9) 



In [12] the relations between the ALH and the Davey-Stewartson equation (together with 
the Ishimori model) have been derived. One can find there expressions for the corresponding 
finite genus solutions. 

The last example stems from the fact that for any fixed n the quantity 



u = r n _ip n q n+ i 
solves the Kadomtsev-Petviashvili (KP) equation, 



(5.10) 



0i (4<9 3 w + dmu + \2udiu) = 3d 22 u (5.11) 
Thus, the results of the previous section yield the following finite genus solution for the KP: 



6 (Ckp( z u z 2, z a) -V^j9 (Ckp( z i, z 2, h) + V 

° 2 (Ckp( Z 1: Z 2,Z 3 )) 



(5.12) 



Ckp( z u z 2 , z 3 ) = C {1) zi + C {2) z 2 + C {3) z 3 + constant (5.13) 

Here I set n = in (|5.10|) and omitted the ^-dependence for all m's as well as the dependence 
on z m for m > 3. This solution differs from the already known one which corresponds to 
an odd-order polynomial V 2g+ i and, what is crucial, has been obtained by considering flows 
near the infinity (ir(P) = oo) which in this case is a ramification (Weierstrass) point ( the 
point £ = oo has only one preimage on the Riemann surface). I cannot at present discuss 
solution fl5.12| ) in detail. For example I do not know whether it is possible to obtain from 
Q5.12 ) any non-trivial real solutions. In any case, solution ( 5.12 ) seems to be interesting and 
worth following studies. 
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To conclude, I want to point out the main differences between the approach of this 
paper and ones used earlier |, |, |, 1. In the IST-based methods the hyperelliptic 
curves appear in the analysis of the spectral data of the scattering problem ( |2.1| ), while the 
dependence on the coordinates is derived from the system (|2.2| ). Here we didn't use the zero- 
curvature representation ( [2.1| )- ( |2.2| ) explicitly (though it is surely hidden in the functional 
equations (|2.9| )- (|2.12| )). We started with an almost arbitrary polynomial (the fact 

that it is related to the transfer matrix of the scattering problem ( [2.1|) wasn't crucial for our 
consideration) and obtained the Zj-, ^-dependence directly from the equations of the ALH 
(and not from the corresponding equations for the transfer matrix T n ). 

As to the algebro-geometrical method, the main distinguishing point is that the approach 
of the work || (and analogous works devoted to other integrable equations) is, so to say 
'global', while ours is 'local'. The authors of used the Baker- Akhiezer function and 
other structures defined for the whole Riemann surface X. At the same time we didn't use 
globally defined objects: each time we introduced some functions depending on the point 
P of the Riemann surface X it was understood that it is defined in some vicinity of some 
distinguished point (B or D). We even didn't discussed the question whether our functions, 
say p QR (P), are well-defined, or single-valued, for all P's, P G X. All we needed is the 
Taylor expansions, say (|4.28|), hence for our purposes it was enough that our functions exist 



locally, for P belonging to some (arbitrary small) neighborhood of the point B (or D). 

At last, I would like to note that the idea to apply the Fay's identity to differential equa- 
tions is far from new. For example, in the book || one can find few examples of how to 
demonstrate that ^-functions of some arguments solve the KP, KdV, sine-Gordon equations. 
However in these examples the main question is how some combinations of differential oper- 
ators (flows) act on a ^-function. To my knowledge the problem of action of these operators 
taken separately hasn't been considered before. Now we know a partial answer: the flows 
near a regular (not Weierstrass) point of a Riemann surface can be described by means of 
the equations of the ALH. Such an appearance of the ALH seems to be new and rather 
interesting. Combined with the results of the works || [K], [IT], [12, |IJ| it can be viewed as 
one more point indicating the 'universality' of this hierarchy. 
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